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Abstract — This paper considers the problem of robust sta- 
bility for a class of uncertain quantum systems subject to 
unknown perturbations in the system coupling operator. A 
general stability result is given for a class of perturbations 
to the system coupling operator. Then, the special case of a 
nominal linear quantum system is considered with non-linear 
perturbations to the system coupling operator. In this case, a 
robust stability condition is given in terms of a scaled strict 
bounded real condition. 

I. Introduction 

An important concept in modern control theory is the 
notion of robust or absolute stability for uncertain nonlinear 
systems in the form of a Lur'e system with an uncertain 
nonlinear block which satisfies a sector bound condition; 
e.g., see [1]. This enables a frequency domain condition for 
robust stability to be given. This characterization of robust 
stability enables robust feedback controller synthesis to be 
carried out using H°° control theory; e.g., see [2]. In a 
recent paper [3], classical results on robust stability were 
extended to the case of nonlinear quantum systems with non- 
quadratic perturbations to the system Hamiltonian. The aim 
of this paper is to extend classical results on robust stability 
to the case of nonlinear quantum systems with nonlinear 
perturbations to the system coupling operator. 

In recent years, a number of papers have considered the 
feedback control of systems whose dynamics are governed 
by the laws of quantum mechanics rather than classical 
mechanics; e.g., see [4]-[16]. In particular, the papers [13], 
[17] consider a framework of quantum systems defined in 
terms of a triple (S,L,H) where S is a scattering matrix, 
L is a vector of coupling operators and H is a Hamiltonian 
operator. The paper [17] then introduces notions of dissipa- 
tivity and stability for this class of quantum systems. In this 
paper, we build on the results of [17] to obtain robust stability 
results for uncertain quantum systems in which the quantum 
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system coupling operator is decomposed as L = L\ + L 2 
where L\ is a known nominal coupling operator and L 2 is 
a perturbation coupling operator, which is contained in a 
specified set of coupling operators W. 

For this general class of uncertain quantum systems, a 
general stability result is obtained. The paper then considers 
the case in which the nominal quantum system (S,Li,H) 
is a linear quantum system in which the Hamiltonian H is a 
quadratic function of annihilation and creation operators and 
the coupling operator L\ is a linear function of annihilation 
and creation operators; e.g., see [7], [8], [10], [11], [16]. In 
this special case, a robust stability result is obtained in terms 
of a scaled frequency domain condition. 

The remainder of the paper proceeds as follows. In Sec- 
tion EH we define the general class of uncertain quantum 
systems under consideration. In Section |IIIJ we consider a 
special class of non-linear perturbation coupling operators. 
In Section ITVl we specialize to the case of a linear nominal 
quantum systems and obtain a robust stability result for this 
case in which the stability condition is given in terms of 
a strict bounded real condition dependent on three scaling 
parameters. In Section [Vj we present some conclusions. 

II. Quantum Systems 

We consider open quantum systems defined by parameters 
(S,L,H) where L = L± + L 2 ; e.g., see [13], [17]. The 
corresponding generator for this quantum system is given by 



G(X) = -i[X,H] + C L (X) 



(1) 



where C L {X) = \L*[X,L) + \[L*,X]L. Here, [X,H] = 
XH — HX denotes the commutator between two operators 
and the notation * denotes the adjoint of an operator. Also, 
H is a self-adjoint operator on the underlying Hilbert space 
referred to as the system Hamiltonian. L\ is the nominal 
system coupling operator and L 2 is referred to as the 
perturbation coupling operator. Also, S is a unitary matrix 
referred to as the scattering matrix. Throughout this paper, 
we will assume that S = I. The triple (S, L, H), along with 
the corresponding generator define the Heisenberg evolution 
X(t) of an operator X according to a quantum stochastic 
differential equation 

dX = (£ L (X)-i[X,H])dt. 

+dA*S^[X,L] + [L*,X]SdA 
+tr [(&XS - X) dk T ] ; 

e.g., see [17]. Here, in the case of operators, the notation 
t denotes the adjoint transpose of a vector or matrix of 



operators; see also [17] for a definition of the quantities dA, 
dA* , and dA which will not be further considered in this 
paper. Also, in the case of standard matrices, the notation < 
refers to the complex conjugate transpose of a matrix. 

The problem under consideration involves establishing 
robust stability properties for an uncertain open quantum 
system for the case in which the perturbation coupling 
operator is contained in a given set W. The main robust 
stability results presented in this paper will build on the 
following result from [17]. 

Lemma 1 (See Lemma 3.4 of [17].): Consider an open 
quantum system defined by (S, L, H) and suppose there ex- 
ists a non-negative self-adjoint operator V on the underlying 
Hilbert space such that 



a V £ V and real constants c > 0, A > 0, t\ > 0, r 2 > 



G(V) + cV < X 



(2) 



where c > and A are real numbers. Then for any plant 
state, we have 

(V(t)) <e- ct (V) + ~, Vi >0. 

c 

Here V(t) denotes the Heisenberg evolution of the operator 
V and (•) denotes quantum expectation; e.g., see [17]. 

A. Commutator Decomposition 

Given a set of non-negative self-adjoint operators V and 
real parameters 7 > 0, 8\ > 0, 82 > 0, £3 > 0, we now 
define a particular set of perturbation coupling operators 
W\. This set Wi is defined in terms of the commutator 
decomposition 

[V,L 2 ] = WF,C]-^2[C,M 0) 

for V £ V where w\, w 2 and £ are given scalar operators. 
We say L2 £ Wi if the following sector bound condition 
holds: 

L* 2 L 2 < 4, C*C + Si (4) 

T 



and 



w\w\ < 82, 
W2W2 < 83. 



(5) 



(6) 



Here, we use the convention that for operator inequalities, 
terms consisting of real constants are interpreted as that 
constant multiplying the identity operator. 
Then, we define 



L 2 : 3wi , W2 , C sucn that ©, (0 and © 
are satisfied and ® is satisfied W £ V 



(7) 

Using this definition, we obtain the following theorem. 

Theorem 1: Consider a set of non-negative self-adjoint 
operators V and an open quantum system (S, L, H) where 
L = L\ + L2 and L2 £ Wi defined in (0. If there exists 



0. 



, T5 > such that \x = — \ [C, [V, Q] is a constant and 



-i[V,H} + C Ll (V)+(^- + 
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, [v,cmc] 



V 2 7 2 2 7 2 
+cV < A, 



then 

(V(t)) < e- ct {V) + 



2- 



L\L X 



[V,£i]*[y,£i] 
2r| 



2r r - 



fJL li- 



ft) 



8l 



for all t > 0. 

Proof: Let V £ V be given and consider Q(V) defined in 
([TJ. Then using (O and the fact that V is self-adjoint, 



-i[V, H] + - (Li + L- 2 )* [V, Li + La] 



1 



[(Li + L 2 )* , V] (L1+L2) 



2 

-i[V,H] 

+ ±(L* 1+ L* 2 ) ([V, L 1 }+w 1 [V, C] + liw 2 ) 
+ \ ([V, LJ* + [V, CYwt + ii*w* 2 ) (Li + L 2 ) 

+iL*[y J L 1 ] + i[v,LirL 2 

+±L* 2 w 1 [V,C} + \[V > t}*w* 1 L 2 
+^L 2 w 2 + ^/j,*w 2 L 2 . 



(9) 



Now 







T1L1 



7"! 



[V,C}*w* lWl [V,C) 



and hence 



+ 



L^^ci + i^cr^i 

2 T *T , [V,(]*W*lWl[V,(] 



< T 1 L 1 Li-\ 2 



(10) 



using ©. Also, 



< (r 2 Ll - ^i) (r 2il - ^ 
= t 2 L\L\ — [iL\w 2 — [i*w 2 L\ 

fl*[lW 2 W 2 



and hence 



[iL 1 w 2 + [J, w 2 Li < ^L^i^ 2 

T 2 

< T 2 LjLl H — 



using (O. Also, 



T"3 / \ T 3 

r 3 2 L*i 2 - L*[V,Lx] - [V,ii]*i 2 
, [V,Li]*[V,Li] 



(ID 



using (01 and ©. 

Substituting (0, (Qj), <H3, O and (O into (0, it 
follows that 



Q{V) < -i[V,H]+C Ll {V) + -l-LtL 1 



S 2 [V,Q*[VX] 



2 1 



2r| 



2r| 



2r| 



-*1 



2 7 2 ^ * ' 2" 1 ' 2r| 



g(V)+cV < A+( ^ + ^ ) + + 



2Tf 



(15) 



Then it follows from (0 that 

r 2 r 2 t 2 " 

'2 ^'5 / y 2 2 2 

Then the result of the theorem follows from Lemma Q] □ 



and hence 



L£[V,Li] + [V;LirL 2 



< T 3 L 2 L 2 -\ j 



^ r2^C + T 3 <5 1 + - 2 (12) 

7 T 3 



using ©. Also, 

< r 4 L 2 ] ( t 4 L 2 



T4 / \ t 4 
r 4 2 L^ 2 - L>![y, C] - [V, C]*<ia 



and hence 



Ll Wl [V,Q + [V,Q*wlL 2 

< T A L 2 L 2 H 2^ 



7 2 



WOW] 



using © and (01. Also, 



= rf L 2 L 2 - iiL* 2 w 2 — fi*w 2 L 2 

fJ,*flW 2 W 2 



and hence 



fj,L 2 w 2 + fi*w 2 L 2 
< T 5 L 2 L 2 ^r 



T f 2 CC + rlS^ 



(13) 



(14) 



III. Non-linear Perturbation Coupling 
Operators 

In this section, we define a set of non-linear perturbation 
coupling operators denoted W 2 . For a given set of non- 
negative self-adjoint operators V and real parameters 7 > 0, 
5i > 0, 6 2 > 0, S3 > 0, consider perturbation coupling 
operators defined in terms of the following power series 
(which is assumed to converge in some suitable sense) 



l 2 = /(c) = J2 Sk < k 

k=Q 



k=0 



S k L k . 



(16) 



Here £ is a scalar operator on the underlying Hilbert space 
and L k = Q k . 
Also, we let 



00 



k=l 



»k-l 



/"(c) = E fc ( fc - l ) s ^ k ~ 

k=l 

and consider the sector bound condition 

/(C)*/(C) < ^CC + Si 

7 

and the conditions 

/'(C)*/'(C) < h, 
/"(C)*/"(C) < s a . 

Then we define the set W 2 as follows: 
L 2 of the form ( TToT ) such that 



(17) 
(18) 

(19) 

(20) 
(21) 



conditions ( fl9] l, ( f20b and (|2TT > are satisfied 

(22) 

In this section, the set of non-negative self-adjoint opera- 
tors V will be assumed to satisfy the following assumption: 



Assumption 1: Given any VeP, the quantity 



where M € C 



2nx2n 



is a Hermitian matrix of the form 



m = ~[c, [v;c]] = -ic[v,c] + |[v,c]c 

is a constant. 

Lemma 2: Suppose the set of self-adjoint operators V 
satisfies Assumption Q] Then 

W 2 c Wi. 

Proof: First, we note that given any V E V and fc > 1, 

^C = [^,C] + C^; 

y( fe = ^c n_1 [^C]C^ n + C fe ^- (23) 

n=l 

Also for any n > 1 such that n < fc, 

[V,C]C = C[^,C]+2m; 

[y, (]( fe -« = ( k - n [V, C] + 2(fc - n)C fc ~""V (24) 

Therefore using d23l l and d24l >. it follows that 

fc 

i/c fc = 5Ic n ~ 1 c A:_n [v;c] + 2(fc-7i)c n ~ 1 c fc ~™~V 

n=l 

+c fe ^ 

fc 

= £<* _1 [v;C] + 2(fc-n)c fc "V + C fc ^ 



= fcc [v, C] + - i)C a* + C n 

which holds for any fc > 0. 

Now given any L 2 £ W 2 , fc > we have 

[v,L k ] = fcc^tv.ci + ^fc-iX^V 

Therefore, 

00 

[V,L 2 ] = X)5*[^ifc] 



(25) 



fc=0 



Now letting 



= /W,C] + /"(C)m- 

Wl = /'(C), and W2 = /"(C), 



(26) 



(27) 



it follows that condition (O is satisfied. Furthermore, condi- 
tions ©, ©, © follow from conditions (O, (|20|), (f2Tb 
respectively. Hence, L 2 £ Wi. Since, L 2 G W 2 was 
arbitrary, we must have W 2 C Wi. □ 

IV. The Case of a Linear Nominal System 

We now consider the case in which the nominal quantum 
system corresponds to a linear quantum system; e.g., see [7], 
[8], [10], [11], [16]. In this case, we assume that H is of the 
form 



H = - [ at a T 1 M 
2 L J 



,# 



(28) 



M 



Mi M 2 



M 2 # 



M 



it 



and Mi = Ml, M 2 = 
annihilation operators on the underlying Hilbert space and 
a# is the corresponding vector of creation operators. In the 
case vectors of operators, the notation * refers to the vector 
of adjoint operators and in the case of complex matrices, this 
notation refers to the complex conjugate matrix. 

The annihilation and creation operators are assumed to 
satisfy the canonical commutation relations: 

it 



M2 ■ Here a is a vector of 



a 






a 


a* _ 




. a * 




a 




# 


( 


. a * . 





where J 



I 
-I 

We also assume L\ is of the form 



= J 

; e.g., see [9], [14], [16]. 



(29) 



L 1 =[N 1 N 2 



a 




a 


. a * . 


= N 


. a * . 



(30) 



where Ni £ 


(Ql xn 


and N 2 £ C 




' L x ' 


= N 


a 













Also, we write 



Ni N 2 



N* NI 



it 



In addition we assume that V is of the form 



a 



v=[ 



a 



(31) 



where P G 
of the form 



is a positive-definite Hermitian matrix 



P 



Pa Pi 



P* 



P 



(32) 



Hence, we consider the set of non-negative self-adjoint 
operators "Pi defined as 



Vi = 



V of the form ( t3Tb such that P > is a 
Hermitian matrix of the form d32l ) 



(33) 

In the linear case, we will consider a specific notion of 
robust mean square stability. 

Definition 1: An uncertain open quantum system defined 
by [S, L, H) where L = L\ + L 2 with L\ of the form d30l l. 
L 2 € W, W is any given set, and H of the form (f28b is said 
to be robustly mean square stable if there exist constants 
Ci > 0, c 2 > and C3 > such that for any L 2 G W 





a(t) 


t 


a(t) 


( 






a#(t) 



< 



a 

a* 



t r 



a 



c 3 Vt > 0. (34) 



Here 



o(t) 
a*(t) 



denotes the Heisenberg evolution of the 



vector of operators 



a 

n* 



; e.g., see [17]. 



In order to address the issue of robust mean square 
stability for the uncertain linear quantum systems under 
consideration, we first require some algebraic identities. 

Lemma 3: Given V G V\, H defined as in d28l ) and L\ 
defined as in d30b , then 



[V,H] = 



[at a T }P 



a 

a* 



a 



a 

a* 



[PJM - MJP] 



a 

a* 



Also, 



\l\[V,L 1 ] + \[lIv]L 1 



= Tr ( PJN* 



1 

~2 

Furthermore, 



a 



{nKJNJP + PJN f JN) 



a 

a* 



[ 



,t 



a 



2JP 



a 



Proof: The proof of these identities follows via straightfor- 
ward but tedious calculations using d29l ). □ 
We now specialize the results of Section [TT] to the case of 
a linear nominal system with L2 G W2 where W2 is defined 
as in Section [HI] In this case, we define 

E 2 a* 



c 



[Eh. E 2 ] 



a 




a 


. a * 


= E 


. a * . 



(35) 



where ( is assumed to be a scalar operator. Then, we show 
that a sufficient condition for robust mean square stability 
is the existence of constants n > 0, T3 > 0, and T4 > 
such that the following scaled strict bounded real condition 
is satisfied: 
1) The matrix 



F 



-iJM - -JN^JN is Hurwitz; 
2 



2) 



G{sI-Fy l B 



< 1 



where 



C 



7 



E 



and 



B = 



T3 



(36) 



(37) 



(38) 



(39) 



This leads to the following theorem. 



Theorem 2: Consider an uncertain open quantum system 
defined by (S, L, H) such that L = L\ + L2 where L\ is 
of the form (|30]l, H is of the form (|28]l and L 2 G W 2 . 
Furthermore, assume that there exists constants T\ > 0, T3 > 
0, and T4 > such that the strict bounded real condition 
(l36l l. (f37T > is satisfied. Then the uncertain quantum system is 
robustly mean square stable. 

In order to prove this theorem, we require the following 
lemma. 

Lemma 4: Given any V G V\, then 



/' 



~[C. \vA\ = -\et.jpje t . 



Hence, the set 



which is a constant. Here, £ = 

j 

of operators V\ satisfies Assumption Q] 
Proof: The proof of this result follows via a straightforward 
but tedious calculation using d29l . □ 
Proof of Theorem [2] If the conditions of the theorem are 
satisfied, then d37l i implies 



C_ 
71 



(sI-F) 1 V2B 



< 1. 



Hence, it follows from the strict bounded real lemma that 
the matrix inequality 



F^P + PF 



2PBB^P+ < 0. 



will have a solution P > of the form d32l ; e.g., see [2], 
[11]. This matrix P defines a corresponding operator V G V\ 
as in (l3Tl i. Then using d38l and 
write. 



it follows that we can 



F 1 P 



PF 
P {26 2 



JE^EJ+^JN^Nj) 



P 



-E^E+^N^N < 0. 



2 7 2 

Hence, we can choose t 2 > and T5 > sufficiently small 
so that 



i?tp. 



PF 
P (262 



-4 + -4 ] J^SJ + ^JN^NJ ] P 



1 '4 / 





t 




a 

. a * . 




a 



2 7 2 

Now, it follows from 



c*c = 

Also, it follows from Lemma [3] that 

[y,C] = -2^JP 

Hence, 



-^N^N < 0. 



that we can write 



(40) 



(41) 



PJ&EJP 



a 



(42) 



Similarly 
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In addition, 

L\Li 



a 



PJN^NJP 



a 

n* 



(43) 



(44) 



Hence using Lemma [3] we obtain 



4[V,H] + C Ll (V)+[ -f 



L\L X 



2t? 
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a* 



2r| 



2 7 2 



2r| 



ftp + pp \ 
1 1 1 ^ PJ&EJP 
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r a +T4_+T B £.t_g 
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ATJ 
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where F 



-iJM 



\JNUN. 

From this, it follows using d40b that there exists a constant 
c > such that condition ([8]) is satisfied with 



A = Tr ( PJN^ 



NJ) > 0. 
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Hence, it follows from Lemma [4] Lemma [2] Theorem Q] and 
P > that 



a(t) 
a#(t) 



t r 
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where A = A 



5 3 

2rr 



Hence, the condition (134-b is satisfied with c\ 
c 2 = c> and c 3 = ^~[F] - °- 

V. Conclusions 



> o, 



2 



XminlP] 



□ 



In this paper, we have considered the problem of robust 
stability for uncertain quantum systems with non-linear per- 
turbations to the system coupling operator. The final stability 
result obtained is expressed in terms of a strict bounded real 
condition. Future research will be directed towards analyzing 
the stability of specific nonlinear quantum systems using the 
given robust stability result. 



